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EPR studies of a Ni4 single-molecule magnet yield the zero-field-splitting (zfs) parameters, D, B1 
and _B|, based on a giant spin approximation (GSA) with 5 = 4. Experiments on an isostructural 
Ni-doped Zn4 crystal establish the Ni'' ion zfs parameters. The 4 
arise from the interplay between the Heisenberg interaction, Jsi • S2, and the 2 
anisotropy, giving rise to mixing of higher lying S 7^ 4 states into the 5 = 4 state. Consequently, J 
directly influences the zfs in the ground state, enabling its direct determination by EPR. 
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PACS numbers: 75.50.Xx, 75.60.Jk, 75.75.+a, 76.30.-v 



The [Ni(hmp)(i?0H)Cl]4 molecule (abbreviated NU) 
possessing the i?OH = dmb ligand {mf^^ 0,11111113) 
represents a model system for carefully examining the 
validity of the giant spin approximation (GSA) which 
has been widely applied in the study of single-molecule 
magnets (SMMs) The GSA assumes the total spin, 
S, of the molecule to be a good quantum number, and 
then models the lowest-lying (25-1-1) magnetic sub-levels 
in terms of an effective spin Hamiltonian of the form: 



H = DSl + Biol 
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The first three terms parameterize anisotropic magnetic 
interactions which lead to zero-field-splitting (zfs) of the 
ground-state multiplet (see red lines in Fig. 1 for the case 
of 5* 4), e.g. spin-orbit coupling, dipolar interactions, 
etc; here, we consider only 2"^^ and 4"^-order operators 
(see for definitions) which are compatible with the S'4 
symmetry of the Ni4™'^ SMM. The final term represents 
the Zeeman interaction associated with the application 
of a magnetic field, B, where g is the Lande g-tensor. 

SMMs are defined by a dominant 2"'^-order uniaxial 
anisotropy, DS^, with Z? < (Qj. Nevertheless, weaker 
4*'^-order terms have been shown to play a crucial role in 
the quantum dynamics of several high-symmetry SMMs 
(especially Mni2-acetate) H H 0, even though the 
precise origin of these terms has not previously been un- 
derstood [ill. In this letter, we show that higher order 
terms [0(2n), n > 1] arise naturally in the GSA through 
the interplay between intrinsic magneto-anisotropy (at 
the sites of individual magnetic ions in the molecule) and 
inter-spin-state mixing (controlled by exchange). These 
findings raise questions concerning the validity of the 
GSA, particularly in terms of its predictive powers. 

The mf"^ SMM is particularly attractive for this in- 
vestigation. The four s = 1 Ni^^ ions reside on opposing 
corners of a shghtly distorted cube (Fig. 1 inset) U, |^ . 
DC susceptibility data (xmT) indicate a relatively large 
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FIG. 1: Field dependence of the 81 eigenvalues corresponding 
to the four-spin Hamiltonian (Eq.|5J. The lowest 9 levels (red 
lines) can be modeled by the GSA with S = 4 (Eq.[TJ. The 
inset shows a schematic of the cubic core of the Kif^^ SMM 
(the blue spheres represent O). 



ground state spin of S* = 4 for the molecule, and a reason- 
able separation (~ 35 K) between this and higher lying 
states with S* < 4 |^, 0] . These properties can be ratio- 
nalized in terms of pure ferromagnetic coupling between 
the Ni^^ ions. In addition, efforts to fit low-temperature 
electron paramagnetic resonance (EPR) and magnetiza- 
tion data to the GSA (Eq Jlland red lines in Fig. ^ have 
been highly successful fllili. Thus, Ni^""'^ displays all 
of the hallmarks of a SMM, yet it exhibits unusually fast 
magnetic quantum tunneling (MQT) at zero field • 



The GSA assumes S to be rigid, thereby ignoring 
the internal magnetic degrees of freedom within a SMM 
which can give rise to couplings to higher-lying states 
(5*— mixing [ill. H^) that may ultimately influence MQT. 
A more physical model, which takes into account zfs 
interactions at the individual Ni^^ sites, as well as the 
exchange coupling between individual magnetic ions, is 
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given by the following Hamiltonian 
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Here, the Sai {a = x, y, z) represent spin projection oper- 
ators, and Si the total spin operator, corresponding to the 
individual Ni'^ ions; di (< 0) and e, respectively param- 
eterize the uniaxial and rhombic zfs interactions in the 
local coordinate frame of each Ni^^ ion; likewise, g j rep- 
resents the Lande g-tensor at each site; finally, the Jy 
parameterize the isotropic exchange couplings between 
pairs of Ni^^ ions. 

For Ni4™'' , the dimension of the four-spin Hamiltonian 
matrix (Eq.|2) is just [(2s-hl)'']^ = 81x81, which is easily 
handled by any modern PC (in contrast to Mni2-acetate 
which has dimension ~ 10^ x 10® More impor- 

tantly, the 3x3 Hamiltonian matrix associated with a 
single Ni^^ ion contains only two zfs parameters, di and 
(in addition to ^J. Furthermore, due to the high sym- 
metry of the molecule, these matrices are related sim- 
ply by the ^4 symmetry operation, and the number of 
exchange constants reduces to just two (Ji and J2, see 
Fig. n inset). Consequently, the four-spin model con- 
tains only a hand full of parameters, each of which can 
be determined independently, often by more than one 
method 0,0,0,1111. Fig. [T] displays the 81 Zeeman-spht 
eigenvalues corresponding to the four-spin Hamiltonian 
(Eq. I2J), using parameters obtained from fits described 
later. The lowest nine levels are fairly well separated 
from higher lying states; these levels, which dominate 
the EPR spectrum, can be equally well described by the 
Hamiltonian corresponding to Eq.^with S = A 
Roughly 20 cm~^ above this ground state multiplet exists 
a grouping of 21 levels which can reasonably be treated 
as three separate 5 = 3 multiplets. There is then a gap 
to a more-or-less continuum of levels. The notion of a 
well defined spin quantum number becomes tenuous at 
this point. 

There are a number of other important reasons why we 
chose to focus on the Ni4™'' member of the Ni4 family. 
To begin with, Ni^^ is readily amenable to substitution 
with non-magnetic Zn. Thus, one can synthesize crystals 
of Znf^^ lightly doped with Ni" The result is a small 
fraction of predominantly s = 1 ZnaNi magnetic species 
diluted into a non-magnetic host crystal. X-ray stud- 
ies indicate that the structures of the Ni4™'^ and Znf™'' 
complexes are virtually identical. Thus, EPR studies of 
the doped crystals provide very reliable estimates of the 
single-ion zfs parameters for Ni^^ in the parent Ni4™'' 
compound (di , and in Eq. I^J . Another remark- 
able feature of the Ni4'°'^ member of the Ni4 family is that 
its structure contains absolutely no solvent of crystalliza- 
tion [^0,H. This is quite rare among SMMs, resulting 
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FIG. 2: (a) ms dependence of the zfs energies [between ms 
and (mg + 1)] within the ground state multiplet. The dashed 
curves show the zfs obtained from Eq. |21 as a function of J. 
The inset defines the Euler angles relating the Ni^' and molec- 
ular coordinates Q. (b) Difference between the data in (a) 
and the J = —333 cm~^ curve, emphasizing the non-linear 
ms dependence of the zfs energies. 



in the removal of a major source of disorder. Indeed, we 
believe that this is the primary reason why the Ni^™*^ 
complex gives particularly sharp EPR spectra 0, 0| . 
In contrast, all of the other solvent containing Ni4 com- 
plexes exhibit rather broad EPR peaks jl| . Details of the 
experimental procedures, including representative EPR 
spectra, are presented elsewhere 0, 0, |l3 • 

We begin by reviewing the results of single-crystal 
high-frequency EPR studies of Kif^^ Q, 0]. Based 
on an analysis using the GSA (Eq. HI), the lowest-lying 
5 = 4 multiplet is split by a dominant axial zfs inter- 
action with D = —0.589(2) cm^^. In the absence of 
higher-order terms, this interaction produces a quadratic 
dependence of the {2S + 1) zero-field eigenvalues on the 
quantum number ms, representing the projection of the 
total spin onto the easy-axis of the molecule. Conse- 
quently, the zfs between successive ms levels should be 
linear in ms- It is these splittings that one measures 
in an EPR experiment, albeit in a finite magnetic field. 
However, using a multi-frequency approach, one can ex- 
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FIG. 3: Angle-dependence of the splitting of the lowest en- 
ergy doublet (mg = ±4 in zero-field) as a function of the field 
orientation within the hard plane. The thin curves are simu- 
lations using Eq.|5|with different J values (see Fig.|5|for color 
codes). The black and red data are the best fits to Eqns. Q 
and|51 respectively. 



trapolate easy-axis data {B \\ z) to zero-field, yielding ac- 
curate determinations of these splittings ; these energy 
spacings are plotted versus in Fig. |21for Ni4'"''. As 
can be seen, the dependence of the zfs values on ms is not 
linear. One can obtain agreement to within experimen- 
tal error by including the 4"^-order axial zfs interaction 
B^Ol (oc S^) in the GSA, with = -1.2 x lO"'^ cm"! 
(black data in Fig. [21 Ij). The S'^ operator produces quar- 
tic ms corrections to the zero-field eigenvalues and, thus, 
cubic corrections to the zfs, as seen in the figure. Unlike 
the 2"'^-order term, which can easily be understood as 
originating from the 2"'^-order zfs interactions at the in- 
dividual Ni^^ sites, the 4*'^-order term in the GSA does 
not have any obvious physical meaning. In this sense, 
it is nothing more than an adjustable parameter in an 
effective model (Eq. As we will see below, this non- 
linear ms dependence of the zfs values is directly related 
to S'-mixing jl ij . 

The four-fold (84) symmetry of the Ni^'"'^ molecule 
forbids 2"'^-order zfs interactions which break axial sym- 
metry. Indeed, we find no evidence for such interac- 
tions based on EPR experiments conducted as a func- 
tion of the field orientation within the hard plane. How- 
ever, a very pronounced four-fold modulation of the spec- 
trum is observed, which can be explained by the 4**^- 
order BfOj [= ^Bi{S^ + S^)] term in the GSA, with 
Sl = ±4 X 10"'' cm~^ 0. Although this interaction is 
allowed by symmetry, it is not obvious how it relates to 
the underlying anisotropy associated with the individual 
Ni^^ ions. Nevertheless, it does explain the fast MQT 
observed in this and other Ni4 complexes |3, 0] . When 



treated as a perturbation to the axial zfs Hamiltonian, 
{S^ + Sy) connects states that differ in ms by ±4 in first 
order and, therefore, lifts the degeneracy of the lowest 
lying ms — ±4 states in 2nd order, leading to a tunnel 
splitting of order 10 MHz. This is an extremely large 
intrinsic tunnel splitting in comparison to other SMMs, 
and can be understood as arising because of the coinci- 
dence of the multiplicity of the ground state (25' -|- 1 = 9) 
and the four-fold symmetry, which gives rise to a leading- 
order off-diagonal zfs interaction which is fourth order in 
the spin operators, i.e. {S'^ -\- Sy) is extremely effective 
at connecting the ms — ±4 states. 

We now attempt to understand the physical basis for 
the existence of the axial and transverse 4*''-order zfs in- 
teractions (-B4O4 and BfOf) deduced on the basis of the 
GSA. From previous studies of a Ni doped Zn4™'' crys- 
tal, we determined not only the zfs parameters associated 
with the Ni^^ ions, but also the orientations of the local 
magnetic axes associated with these interactions relative 
to the crystallographic axes 0. However, the key point 
is that the Hamiltonian matrices for the individual s = 1 
Ni^^ ions have dimensions 3x3. Therefore, terms ex- 
ceeding 2"'^-order in the single-spin operators (sf^, s^y, 
etc.) are completely unphysical. If one assumes that 
the ground state for the Ni4™'' molecule corresponds to 
a rigid 3 = 4 spin, one can then project the single-ion 
anisotropics onto the 5 = 4 state using irreducible ten- 
sor operator methods 0. However, after rotating from 
local to molecular coordinates, the projection is nothing 
more than a summation of the individual zfs matrices. 
Consequently, such a procedure does not produce terms 
of order four in the spin operators 0. Therefore, the 
need to include 4"^-order zfs interactions in an analysis 
of the EPR data for Ni4™'' may be taken as evidence 
for a breakdown of the GSA. We note that agreement 
in terms of the 2"'^-order parameters is very good. In 
particular, the molecular D value agrees to within 10% 
with the value obtained from projection of the single-ion 
anisotropics onto the 5 = 4 state In addition, al- 
though the single-ions experience a significant rhombic 
zfs interaction (e/d ^ 0.23), symmetry considerations 
guarantee its cancelation when projected onto the 5 = 4 
state. Therefore, this approach is completely unable to 
account for the MQT in Wif'^'°. 

In view of the above, one is forced to use a more re- 
alistic Hamiltonian (Eq. 2) which takes into account all 
spin states of the molecule. The isotropic exchange in- 
teraction, Jij, in Eq. 2 connects states having the same 
spin-projection 12]. Consequently, it does not operate 
between states within a given spin multiplet, it simply 
lifts degeneracies between states with different multiplic- 
ity (see Fig. The addition of anisotropic terms to 
Eq. 2 results in zfs within each multiplet which, in turn, 
gives rise to weaker ms-dependent corrections to the ex- 
change splittings. Thus, we see that J directly modifies 
the zfs energies within a given spin multiplet via interac- 



4 



tions (S'-mixing) with nearby excited spin states. In the 
Umit J » d one can expect these corrections to be negh- 
gible. However, in the present case, where J/d ^ 1.3, one 
can expect these corrections to be significant. Further- 
more, since the corrections involve higher order processes 
whereby the underlying anisotropic interactions feedback 
into themselves via exchange coupling to nearby spin- 
multiplets, it is clear that this will generate 'effective' 
interactions that are 4*''-order (i.e. 2"'^-order squared) in 
the spin operators (as well as higher order terms). How- 
ever, these 4"^-order interactions have no real physical 
basis other than that they arise due to the competing 
isotropic and anisotropic interactions in Eq. 2, resulting 
in S'-mixing. 

The influence of J on the lowest lying (nominally 
5 = 4) multiplet is abundantly apparent in Fig. 2, where 
we compare zfs energies determined via the four-spin 
Hamiltonian (Eq. ^ for different values of the exchange 
interaction strength, with those determined experimen- 
tally (blue data points) and from a fit to the experimen- 
tal data using the GSA (Eq. ^ black data points). The 
magnitudes of d = —4.73 cm~^ and e = —1.19 cm~^ 
were established from combined fits to both easy-axis zfs 
data (red points in Fig. 12), and from hard-plane rotation 
measurements of the four-fold oscillation of the ground 
state splitting (Fig.|3| see also We made one simpli- 
fying assumption by setting Ji = J2 = J, based on DC 
XmT data [l^- Regardless, this in no way invahdates 
the main conclusion of this letter: namely, that J influ- 
ences the ground state zfs through .S'-mixing. The polar 
angle, 9 (see Fig. |21 inset), between the local Ni'^^-ion z- 
axes and the crystallographic z-axis was fixed at 15° on 
the basis of the Ni/Zn studies Q. We additionally in- 
cluded a dipolar coupling (not shown in Eq. between 
the four Ni^^ ions using precise crystallographic data and 
no additional free parameters [llj. The remaining free 
parameters were gx = Qy ~ 2.23, gz = 2.25 and an ad- 
ditional Euler angle (0 = 59°) illustrated in the inset 
to Fig. [3 A more in-depth account of the fitting proce- 
dure will be given elsewhere 0| . The obtained value of d 
agrees to within 12% with the value determined indepen- 
dently from measurements on the Ni-doped Zn4 crystal 
[d ~ —5.30(5) cm~^]; the remaining parameters agree 
to within the experimental error [e = —1.20(2) cm~^, 
.9, = 3^ = 2.20(5), = 2.30(5)]. 

One can clearly see that, by reducing the separation 
between the ground S = 4 multiplet and the lowest ex- 
cited states (by reducing J) , one can reproduce both the 
nonlinear ms dependence of the zfs energies (Fig. 1^ . 
which was attributed to the term in the GSA [l|, 
and the four-fold oscillation of the ground-state splitting 
observed from hard-plane measurements (Fig.|21), which 
was attributed to Bf 2]. This is quite a remarkable re- 
sult, because it implies that one can deduce J directly 
from the spectroscopic information obtained via an EPR 
experiment. Indeed, the value of J = —5.9 cm~^ de- 



termined from these fits is in good agreement with the 
value of —7.05 cm~^ deduced on the basis of fits to xmT 
data to Eq. El 16]. All of the apparent 4*'^-order be- 
havior vanishes if one sets J >> d, as expected in such 
a limit in which the ground state spin value is a good 
quantum number (due to the absence of S'-mixing). In 
the opposite extreme {J ^ —3 cm~^), we start to see 
evidence for even higher order corrections to the zfs en- 
ergies (6"^ order). A cubic polynomial exhibits only one 
turning point (at ms = —0.5 in Fig. EJ, whereas the 
green data in Fig. [3 clearly display more than one turn- 
ing point when one recognizes that all of these curves 
must be antisymmetric about ms = —0.5. Therefore, it 
is apparent that one should not limit the GSA to 4**^- 
order terms for SMMs with relatively low-lying excited 
spin states. In fact, one cannot rule out equally good 
fits to experimental data which include 6*'' and higher- 
order zfs interactions. Consequently, one should be care- 
ful about making predictions on the basis of the GSA, 
particularly at vastly different energy scales compared to 
the experiments used to establish the GSA zfs parame- 
ters (e.g. EPR vs. MQT). Indeed, we find a difference 
of almost a factor of 10 between the ground-state tunnel 
splittings deduced from Eqs. 1 and 2 using the optimum 
zfs parameters for Ni4"*. We note that the situation in 
Ni4 is not dissimilar to many other SMMs, including the 
most widely studied Mni2-acetate, for which similar 4**^- 
order zfs interactions and low-lying excited spin states 
are found [ll[l3- 

Finally, we note that the most unambiguous method 
for estimating exchange couplings in polynuclear metal 
complexes involves determining the exact locations of ex- 
cited spin multiplets. However, the magnetic-dipole se- 
lection rule forbids transitions between states with differ- 
ent multiplicity. Therefore, such an undertaking is usu- 
ally only possible using neutrons 19]. However, Figs. 2 
and 3 clearly show that J can be estimated on the basis 
of zfs of the lowest lying multiplet. Due to the resultant 
S'-mixing, it may be feasible to observe inter-spin-state 
EPR transitions directly via far-infrared techniques. 
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